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ABSTRACT 

We present the result of a search of the Milagro sky map for spatial correlations 
with sources from a subset of the recent Fermi Bright Source List (BSL). The BSL 
consists of the 205 most significant sources detected above 100 MeV by the Fermi Large 
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Area Telescope. We select sources based on their categorization in the BSL, taking all 
confirmed or possible Galactic sources in the field of view of Milagro. Of the 34 Fermi 
sources selected, 14 are observed by Milagro at a significance of 3 standard deviations 
or more. We conduct this search with a new analysis which employs newly-optimized 
gamma-hadron separation and utilizes the full 8-year Milagro dataset. Milagro is sen- 
sitive to gamma rays with energy from 1 to 100 TeV with a peak sensitivity from 10-50 
TeV depending on the source spectrum and declination. These results extend the ob- 
servation of these sources far above the Fermi energy band. With the new analysis 
and additional data, multi-TeV emission is definitively observed associated with the 
Fermi pulsar, J2229. 0-1-6114, in the Boomerang Pulsar Wind Nebula (PWN). Further- 
more, an extended region of multi-TeV emission is associated with the Fermi pulsar, 
J0634.0-hl745, the Geminga pulsar. 

Subject headings: gamma rays: observations — pulsars: general — supernova remnants 



Introduction 



The Milagro gamma-ray observatory has perf ormed the most sen sitive survey of 1 to 100 TeV 

2007al lb[). The Milagro data set is 



gamma rays from the Northern Hemisphere sky (jAbdo et al 



ideal for searching for new classes of g amma-ray sources . The recent release of the Bright Source 
List (BSL) by the Fermi collaboration (jAbdo et al.ll2009l ) presents such an opportunity by looking 
for coincidences of > 1 TeV emission with these GeV sources. There are 34 sources in the BSL 
within Milagro's field of view that are not associated with extragalactic sources. We present a 
search of the Milagro data for excesses between 1 and 100 TeV coincident wi th these 34 potent ial 
Galactic sources. The analysis presented here differs from previous analyses ( Abdo et al. 2007al lbl. 
20081 ) by optimizing the event weighting and Gaussian weighting separately in bins of event size 
(measured with the fraction of channels hit in the instrument). With the improved analysis and 
an additional year and a half of data, the sensitivity has increased by 15% to 25%, depending on 
the spectrum of the source. 



2. Analysis and Results 

We select Fermi-LAT sources in the field of view of Milagro (with 5 > -5°) based on their 
categorization in the BSL. Sources are selected which are confirmed or potential Galactic sources. 
Sources that are identified as extragalactic are omitted. Sixteen of the selected sources were cate- 
gorized in the BSL as confirmed pulsars (PSR) and one is a high-mass X-ray binary (HXB). Five 
sources have a potential association with an SNR, and 12 have no clear association. For each of 
these 34 sources, we calculate the statistical significance of the Milagro data at the BSL position 
and estimate the fiux or fiux limit under the assumption that the emission is from a single point 
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source. 



The flux measurements given in Table [T] are derived with a similar approach to (jAbdo et al 

2007bl ). The flux is measured with an assumed spectrum of E""^'^ without a cutofl'. The dependence 



of the calculated flux on the true spectrum is minimized when the flux is quoted at the median en- 
ergy of the hypothesized spectrum. The median energy depends on Declination and varies between 
32 and 46 TeV for 6 in the range of 10° and 60°. At a Declination of -5°, the median energy of the 
hypothesized spectrum is 90 TeV. We quote the flux for all sources above 3(T at a representative 
value of 35 TeV. It should be noted that the median energy used is for the assumed spectrum and 
not experimentally measured. In particular, a source may in fact cut off before 35 TeV (the Crab 
for example) and our analysis would still report a flux at 35 TeV. The energy spectrum of each 
reported source is the subject of a paper in preparation. 

The results of this search are summarized in Tabled) Of the 34 targets, 14 have a significance 
greater than So". Six o f these are associate d with sources or candidates from the first Milagro survey 
of the Galactic plane \khdo et ahlboOTbl l. The Crab, MGRO 2019+37, MGRO 1908+06, MGRO 
2031+41, and Milagro candidates C3 (likely associated with Geminga) and C4 (likely associated 
with the Boomerang PWN) are all near LAT GeV sources. In the Milagro data set, the 3(T-5(T 
observations are fairly marginal because they cannot be convincingly discerned from background 
when statistical penalties for searching the entire sky are taken into account. However, with LAT 
points as a trigger for the search, the statistical penalties are reduced. The probability of a single 
3(T false-positive in 34 samples of pure background is only ~4.4%. The probability of 4 or more 
excesses at or above 3(7 in 34 trials is ~ 1.5 x 10~^. It is very likely that most of our 3(T excesses are 
due to multi-TeV emissior0. We, therefore, see strong evidence for multi-TeV emission associated 
with Galactic LAT BSL sources as a class, even if individual sources are not strong enough to 
definitively distinguish. 

There is some contribution to these measurements from the Galactic diffuse emission, but that 
contribution is small. We can make a co nservative estimate by taking the Milagro measurement 
of the diffuse emission (lAbdo et al.ll2008l ) at its highest value, in the inner Galaxy (30° < I < 65°, 



\b\ < 2°). Using this value, we expect 5.3 x 10"^^ TeV~^ s cm in a 1° bin at 35 TeV, which is 
only about ~15% contamination for the weakest sources in Table [TJ The GALPROP conventional 
model, for comparison, would only constitute ~3% contamination. The contamination is likely 
lower than suggested by the Milagro measurement becau se of unresolved sources, such as many 
of the sources from Tabled) It has even been suggested (jCasanova &: Dingual2008l ) that most of 
the Milagro diffuse measurement could be due to unresolved sources. Finally, the Fermi points 
observed at 3a in the Milagro data occur near local maxima in the Milagro data. In contrast, the 



^ Alternatively, using the False Discover Rate method ( Miller et al. 2001 : Hopkins et al. 20021 ) and requiring an 
estimate of 1% of the members of the selected candidates to be a false discovery, gives the same list of candidates. 
Changing the contamination fraction criterion from .01 to .001 (or to 0.1), would have included one fewer (or 3 more) 
sources, respectively. 
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diffuse emission is expected to vary slowly across the Galaxy. 



3. Discussion 



From this analysis, it appears quite common for Galactic 100 MeV - 100 GeV sources to have 
associated multi-TeV emission. This association is notable for pulsars, where 9 of 16 pulsars from 
the BSL are on our list of likely multi-TeV emitters. The pulsars in the BSL which have less than 3cr 
significance in Milagro data tend to lie off the Galactic plane. The pulsars off the plane are typically 
older, havi ng traveled far from the ir origin after the kick they received from the initial asymmetric 
supernova (jGunn fc Qstrikeiill970l ). Of the SNR sources on the list, we see 3 of 5. Interestingly, 
we see only 2 of the 12 unidentified sources. These unidentified sources may be extragalactic and 
not visible with this analysis which was optimized for high-energy emission. 

Figure [T] and [2] shows the regions in the Milagro data around the indicated LAT sources. Eight 
of the 13 sources are associated with previously reported >TeV sources or candidates: 

OFGL J0534. 6+2201 is the young Crab Pulsar. Its associated pulsar wind nebula (PWN) 
is a standard reference source in TeV astronomy. 

OFGL J0617. 4+2234 is associated with SNR IC443, which is interacting w ith a nearby 



large molecular cloud. An associated x-ray feature has been interpreted as a PWN (jOlbert et al 



200ll ). implying the existence of a pulsar, bu t the no pulsed ern ission has yet been detected. IC443 
was first reported abov e 1 TeV by MAGIC ([Albert et ahlboO?! ) and later confirmed by VERITAS 
(jHumenskv et al.l |2008| ) . The flux reported in Table [J is somewhat higher than the flux predicted 
by extrapolating the MAGIC fit, but is roughly consistent after allowing for the extremes of the 
statistical and systematic errors of the two measurements. 

OFGL J0634. 0+ 1745 is the Geminga pulsar. Geminga is a rela tively old (342 kyr) but very 
near (169 pc) pulsar (jManchester et al.l l2005l : iHalpern &: Holtl Il992l ). It is the most significant 
Fermi-LAT source in the northern sky, but emission over 1 TeV has only been reported by Milagro 
as candidate C3 with too low a significance to be classified as a definitive detection. Milagro ob- 
serves an emission region that is extended by several degrees as shown in Figure [21 The significance 
reported in Table [1] has been computed assuming point source emission, but if we instead assume 
that the source is due to emission from an extended region and convolve a 1° Gaussian with the 
energy-dependent point spread function, the significance at the location of OFGL J0634. 0+1745 in- 
creases to 6.3(7. The local maximum of the Milagro excess is at RA=6h32m28s, Dec=17°22m. Given 
the high significance, we regard this as a definitive detection of extended emission from Geminga. 
A spatial Gaussian fit to the data yields a region with a standard deviation of 1.30°±0.20°. For 
comparison, the analogous fit for the Crab, which is effectively a point source, has a a of 0.6° This 
suggests that the full width at half maximum of the region of emission in the vicinity of Geminga 
is 2.6^g g°, after accounting for the point spread function. The large extent (implying an emission 
region of some 5 to 10 pc extent) is likely due to the nearness of the source and may arise from 
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a pulsar-driven wind; it is consistent with HESS observation of more distant PWN with an angu- 
lar size of ~10 pc. This may also e xplain why the so urce has not yet been observed by Imaging 
Atmospheric Cherenkov Telescopes (|Maier et al.ll2008l ). 



OFGL J1907.5-K0602 is associated with MGRO .T19n8-^n6 (lAbdo et alJ l2nn7bl ^. This pul- 
sar was discovered by the LAT and is also coincident wi th AGILE source lAGL J1908-I-0613 
(jPittori et alJ bood ) and EGRET source GEV J1907-h557 dLamb fc MacombI Il997l l. The multi- 
TeV emission was first reported by Milagro. HESS both coi ifirmed the Milagro detect ion and was 



also able to identify this source as extended by 0.211^;^^ (jPiannati-Atai et alJ 120071 ') . The peak 



of the Milagro detection occurs at RA=19h6m44s, Dec=5°50m with a 1 sigma error circle of 0.27° 
and a local peak significance of 8.1a. The peak of the Milagro emission is 0.3° from the pulsar, but 
consistent with the pulsar's location within the measurement error. 

OFGL J1923.0-I-1411 is associated with SNR G49. 2-0.7 (W51) which is in a star-forming; 
regio n and near molecular clouds. Recently, a >TeV source, HESS J1923-I-141 (jFeinstein et al 



20091 ). has been detected which is spatially extended and coincident with the Fermi source. 



OFGL J2020.8-I-3649 is associated with MGRO J2019+37. This is the most significant 
source in the Milagro data set apart from the Crab. The young central pulsar has a period of 
104 ms and an estimated age of 17.2 kyr. This source w as also detected by A GILE and EGRET. 
It was AGILE that first identified the GeV pulsations (IHalpern et al.l l2008l ) and that discovery 
was confirmed with Fermi data. The peak of the flux measured by Milagro is at RA=20hl8m43s 
Dec=36°42m with a 0.09° 1-sigma error circle. The position of the excess is ~0.3° from the pulsar. 

OFGL J203 2.2-I-4122 is a LAT ide ntifled pulsar that is spatially coincident with the HEGRA 
source J2032+41 (Aharonian et al.l [20o3). MGRO J20 31+41, and the MA GIC source J2032+4130 
(jAlbert et al.ll2008l ) . The Milagro source was reported (jAbdo et al.ll2007bl ) with an extent of 3° , but 
it appears that the Milagro extended source may be due to two or more overlapping sources with 
a potential additional diffuse contribution from the highly emissive Cygnus region. The location of 
the Milagro peak is RA=20h31m43s and Dec=40°40m with a statistical error of 0.3°. 

OFGL J2229. 0-1-6114 is coincide nt with the radio pulsar J2229-I-6114 which has been pre- 
viously associated dHalpern et ahlboOlh with t he EGRET source 3 EG J2227+6122. The period 



of this pulsar is 52 ms, its distance is 0.8 kpc (IKothes et a 
be 10.5 kyr and E is 2.2 x 10'^'' ergs/sec ( Manchester et al 



200ll'l. and the age is es timated to 
2OO5I : iHalpern et al.l[200ll ). Milagro 



detects a 6.6a excess at the position of the pulsar and a local maximum of 6.8a". The peak of the 
Milagro excess is RA=22h28ml7s Dec=60°29n i with a statistical p osition error of 0.36. This source 
was reported as candidate C4 by Milagro in (jAbdo et al.ll2007bl ). With the additional data and 
improved analysis presented here, this source is elevated to a high-confidence detection. Milagro 
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also identifies this source as clearly not a point source, with a long extension to the soutlj^ 



The remaining five objects with greater than 3a excess in the Milagro data have not been 
previously detected above 1 TeV energies: 



OFGL J0631. 8+1034 is the radio pulsar J0631+10 (jZepka et al.lll996l ). This pulsar has a 
period of 288 ms and an esti mated age of 43.6 kyr, a distance of 6.55 kpc and E of 1.7x10"^^ 
erg/s (iManchester et al.l I2OO5I ) . The VERITAS upper-limit for this region is 1.3% of the Crab 
dMaier et ahlbood ^. 



OFGL J 1844. 1-0335 is unassociated with any known source. It is an interesting source 
because it occurs at a Declination at the edge of Milagro's sensitivity and, if the Milagro observation 
is real, it is extremely bright above 1 TeV. It is in the region of the Galactic plane surveyed by 
HESS (jAharonian et al.ll2006l ) but was not detected. To account for the HESS non-detection, the 
source would have to be extended or have a very hard spectrum extending to high energy. 

OFGL J1900.0-|-0356 has no known associations. 

OFGL J195 4.4-I-2838 is coinci dent with SNR G65. 1+0.6 which has been associated with 
PSR 1957+2831 (|Tian fc Leahvlbood l. 



OFGL J1958 .1+2848 is a LAT-d iscovered pulsar that is associated with the EGRET source 
3EG J1958+2909 (|Hartman et al.lll999l ^. 



OFGL J2021. 5+4026 is a LAT-discovered pulsar that is coincident with the gamma-Cygni 
SNR. This source is located in the Cygnus region that is detected by Milagro as having a broad 
extended excess. 

The relationship between the Fermi and Milagro source fluxes and upper limits for these 34 
sources is shown in Figure [3l The BSL values for the integral flux are shown with the Milagro 
measurements of the differential flux at 35 TeV. The 35 TeV fluxes are roughly correlated with the 
measurements between 100 MeV and 100 GeV but the correlation is not strong, with a correlation 
coefficient of the 3cr points in log space of only 0.2. One possible explanation for the pulsar 
variation is that the pulsed emission is expected to be beai ned (and thus viewing-a ngle dependent) 
and the unpulsed multi-TeV emission is likely unbeamed (jGaensler &: Sland l2006l ) . The spectrum 
that connects the Milagro flux to the Fermi flux is universally softer than 2.0 and closer to 2.3, 
depending on the source. 

We have found that the population of Fermi sources observed at or above 3cr by Milagro is 
dominated by pulsars and/or their associated PWN. Of the 4 high-confidence Milagro detections 



^Note added in press: The Fermi-LAT collaboration has submitted a paper announcing the discovery of a new pul- 
sar - not included in the BSL - with the current best position of RA=339.561, Dec=59.080 (|Fermi-LAT Collaboration 
200' 



^). Milagro observes a 4.7a excess at the location of the pulsar. It may be that the large size of the multi-TeV 



emission associated with OFGL J2229. 0+6114 is in fact due to these two nearby sources. 
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associated with pulsars of known periodicity and distance, 3 (namely J0534. 6+2201, J0634. 0+1745, 
and J2229.0+6114) have E/d^ above 10^^ ergs s~^ kpc~^ where E is the spin down luminosity and 
d is the distance to the pulsar. The distance on th e fourth (J2020.8+3649 ) is uncertain. Using the 
3-4 kpc distance implied by x-ray measurements (jVan Etten et al.ll2008l ) rather than the 12 kpc 
measurement implied by the pulsar dispersion measurement, it too has E/ d^ above 10^^ ergs s ~^ 



kpc . A similar association with high E/d^ pulsars is reported by HESS ( Carrigan et al.l 120071 ). 



Since the pulsed emission is beamed and the PWN is not, all high E /d^ are possible candid ates 
formulti-TeV emission. We have searched the ATNF pulsar database (jManchester et al.ll2005l ) for 
northern-hemisphere pulsars with a high E/d^, which were not reported in the Fermi BSL. Of the 25 
highest E/d^ pulsars, there are 10 in the northern-hemisphere and 5 not identified as GeV sources by 
Fermi. These 5 are J0205+6449, J0659+1414, J1930+1852, J1913+1011 and J1740+1000. Of these, 
the largest statistical significance was 3.3 standard deviations (PSR J1930+1852), not significant 
enough to claim this as new source of multi-TeV gamma rays (though follow-up observations are 
warranted) . 

Finally, it is interesting to n ote that of the sources published in the Milagro survey of the 
Galactic plane (jAbdo et al.ll2007bl ). all 4 sources and two of the 4 source candidates are now strongly 
associated with pulsars, suggesting that most of the Milagro sources are multi-TeV PWN. We also 
note the high efficiency with which MeV to GeV pulsars are observed above 1 TeV, and a qualitative 
picture is emerging where the typical Galactic multi-TeV source is a PWN associated with a MeV 
to GeV pulsar. 
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J06 1 7.4+2234 



JQ63 1.8+ 1034 



JIS44.l^ 




Fig. 1. — The 3a sources from Table [H omitting J0634. 0+1745 (shown in Figm'e[2|) and the Crab. 
Each frame shows a 5°x5° region with the LAT source positions indicated by white dots. The error 
on the Fermi source locations is quite small on this scale, typically between 0.1 and 0.2 degrees, 
depending on the source. The data has been smoothed by a Gaussian of width varying between 
0.4° and 1.0°, depending on the expected angular resolution of events. Horizontal axes show Right- 
Ascension and vertical axes show Declination. The colors indicate the statistical significance in 
standard deviations. 
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J0634.0+1745 J0634.0+1745 




06h40m 06h30m 06h40m 06h30m 



Fig. 2. — The signifcance of the Milagro data in a 5°x5° region around Fermi source J0634. 0+1745 
the Geminga pulsar. The location of the Fermi source is identified by a white dot. The figure 
on the left shows the significance map after smoothing by the Milagro point-spread function. The 
figure on the right shows the same region smoothed by an additional 1° Gaussian in order to search 
for an extended emission region. The color scale shows the statistical significance. 
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Fig. 3. — Flux estimates and upper limits for the 34 Fermi sources. The horizontal axis quotes the 
integral Fermi fiux from 100 MeV to 100 GeV and the vertical axis gives the Milagro fiux or upper 
limit at 35 TeV. Lines are shown with the extrapolation of the Fermi fiux to Milagro energies, 
assuming an i?"^ *^, E~'^-^ and E~'^-^ spectrum. The left panel shows the results for the 16 pulsars 
and the right panel shows the results for the remaining 18 sources. 



Table 1. Summary of the search for TeV emission from sources in the Fermi BSL. 



Name 
(OFGL) 



type 



RA 
(deg) 



DEC 
(deg) 



(deg) 



b 

(deg) 



Flux 



Signif. 

(a's) 



TeV 
assoc. 



J0007.4+7303 


PSR 


1.85 


73.06 


119.69 


10.47 


< 90.4 


2.6 




J0030. 3+0450 


PSR 


7.60 


4.85 


113.11 


-57.62 


< 20.9 


-1.7 




J0240. 3+6113 


HXB 


40.09 


61.23 


135.66 


1.07 


< 26.2 


0.7 


LSI +61 303 


J0357.5+3205 


PSR 


59.39 


32.08 


162.71 


-16.06 


< 16.5 


-0.1 




J0534. 6+2201 


PSR 


83.65 


22.02 


184.56 


-5.76 


162.6 ± 9.4 


17.2 


Crab 


J0613. 9-0202 


PSR 


93.48 


-2.05 


210.47 


-9.27 


< 60.0 


-0.0 




J0617.4+2234 


SNR" 


94.36 


22.57 


189.08 


3.07 


28.8 ± 9.5 


3.0 


IC443 


J0631. 8+1034 


PSR 


97.95 


10.57 


201.30 


0.51 


47.2 ± 12.9 


3.7 




J0633. 5+0634 


PSR 


98.39 


6.58 


205.04 


-0.96 


< 50.2 


1.4 




J0634. 0+1745 


PSR 


98.50 


17.76 


195.16 


4.29 


37.7 ± 10.7 


3.5 


MGRO C3 


















Geminga 


J0643. 2+0858 




100.82 


8.98 


204.01 


2.29 


< 30.5 


0.3 




J1653. 4-0200 




253.35 


-2.01 


16.55 


24.96 


< 51.0 


-0.5 




J1830. 3+0617 




277.58 


6.29 


36.16 


7.54 


< 32.8 


0.2 




J1836. 2+5924 


PSR 


279.06 


59.41 


88.86 


25.00 


< 14.6 


-0.9 




J1844. 1-0335 




281.04 


-3.59 


28.91 


-0.02 


148.4 ± 34.2 


4.3 








282 16 


-1 64 




-0 12 




1 7 




J1855.9+0126 


SNR" 


283.99 


1.44 


34.72 


-0.35 


< 89.5 


2.2 




J1900.0+0356 




285.01 


3.95 


37.42 


-0.11 


70.7 ± 19.5 


3.6 




J1907.5+0602 


PSR 


286.89 


6.03 


40.14 


-0.82 


116.7 ± 15.8 


7.4 


MGRO J1908+06 


















HESS J1908+063 


J1911. 0+0905 


SNR" 


287.76 


9.09 


43.25 


-0.18 


< 41.7 


1.5 




J1923.0+1411 


SNR" 


290.77 


14.19 


49.13 


-0.40 


39.4 ± 11.5 


3.4 


HESS J1923+141 


J1953.2+3249 


PSR 


298.32 


32.82 


68.75 


2.73 


< 17.0 


0.0 




J1954.4+2838 


SNR" 


298.61 


28.65 


65.30 


0.38 


37.1 ± 8.6 


4.3 




J1958.1+2848 


PSR 


299.53 


28.80 


65.85 


-0.23 


34.7 ± 8.6 


4.0 




J2001. 0+4352 




300.27 


43.87 


79.05 


7.12 


< 12.1 


-0.9 




J2020.8+3649 


PSR 


305.22 


36.83 


75.18 


0.13 


108.3 ± 8.7 


12.4 


MGRO J2019+37 


J2021. 5+4026 


PSR 


305.40 


40.44 


78.23 


2.07 


35.8 ± 8.5 


4.2 




J2027.5+3334 




306.88 


33.57 


73.30 


-2.85 


< 16.0 


-0.2 




J2032.2+4122 


PSR 


308.06 


41.38 


80.16 


0.98 


63.3 ± 8.3 


7.6 


TEV 2032+41 


















MGRO J2031+41 


J2055.5+2540 




313.89 


25.67 


70.66 


-12.47 


< 17.6 


-0.0 




J2110.8+4608 




317.70 


46.14 


88.26 


-1.35 


< 24.1 


1.1 




J2214.8+3002 




333.70 


30.05 


86.91 


-21.66 


< 20.7 


0.6 




J2229.0+6114 


PSR 


337.26 


61.24 


106.64 


2.96 


70.9 ± 10.8 


6.6 


MGRO C4 


J2302.9+4443 




345.75 


44.72 


103.44 


-14.00 


< 13.2 


-0.6 





Note. — The source identity in the OFGL catalog is given with the source location in celestial and galactic 
coordinates. We give the measured flux for all sources above 3cr at a characteristic median energy of 35 TeV. The 
2(T upper limits are given for other sources. The statistical significance and nearby TeV associations arc noted. 

^Tho BSL association with a known SNR, is based on similar location. 
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Erratum: Milagro Observations of Multi-TeV Emission from Galactic Sources in the 
Fermi Bright Source List (2009, ApJ, 700L, 127A) 

The position of the peak of the Milagro excess coincident with OFGL J2229. 0+6114 was 

incorrectly reported as RA=22h28ml7s Dec=60°29m. The correct position is RA=22h28m44s 
Dec=61°10m with a statistical position error of 0.165°. 



